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Abstract 

Shell model calculations in the full p/-shell are carried out for the A=50, 51 and 52 
isobars. The most frequently used effective interactions for the p/-shell, KB3 and 
FPD6 are revisited and their behaviour at the N=28 and Z=28 closures examined. 
Cures to their -relatively minor- defaults are proposed, and a new mass dependent 
version called KB3G is released. Energy spectra, electromagnetic transitions and 
moments as well as beta decay properties are computed and compared with the ex- 
periment and with the results of the earlier interactions. A high quality description 
is achieved. Other miscellaneous topics are addressed; the Coulomb energy differ- 
ences of the yrast states of the mirror pair 51 Mn- 51 Fe and the systematics of the 
magnetic moments of the N=28 isotones. 

Key words: A=50, A=51, A=52, Shell Model, Effective interactions, Full p/-shell 
spectroscopy, Level schemes and transition probabilities, Gamow- Teller beta 
decays, Half-lives, Magnetic moments, Coulomb energy differences. 
PACS: 21.10.-k, 27.40.+Z, 21.60.Cs, 23.40.-s 



1 Introduction 



The p/-shell has been the focus of a lot of activity in nuclear structure during 
the last years. Prompted in some cases by the large scale shell model results, 
that indicated the presence of a region of deformation around 48 Cr [1-3], a lot 
of new experiments and calculations have been carried out, addressing many 
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different issues; deformed bands and band termination [4], yrast traps [5], 
high K isomers [6], coexistence of deformed bands of natural and non-natural 
parity [7], effects of neutron proton pairing [8], etc. A very recent highlight 
has been the discovery of an excited deformed band in 56 Ni [9] coexisting with 
the spherical states based in the doubly magic ground state. In addition to the 
exact shell model diagonalizations, the new Monte Carlo techniques, SMMC 
and MCSM have been extensively applied to this region [10-13]. Mean field 
descriptions of various kinds have also been used to explore different issues 
concerning this deformation region [14,15]. 

In this article we extend the full p/-shell calculations up to A=52. Detailed 
results for 50 Cr, 50 Mn and 52 Fe using KB3 have been already published in refs 
[5,16,17] and we will not deal with them here because the new interaction 
KB3G gives equivalent results. In the cases of 51 Cr, 52 Cr, 51 Mn and 52 Mn 
we have carried out the full p/-shell calculation for the yrast states only. To 
perform detailed spectroscopy in the full space would have demanded a huge 
amount of computer time, not justified by the improvement on the results, 
as we have checked. Hence, for the non-yrast states we shall present results 
in a truncated (t=5) space (no more than 5 particles are allowed to excite 
from the lf 7 / 2 subshell). At this truncation level, the most relevant states are 
sufficiently converged. 

As we have discussed in detail elsewhere [3] our usual effective interaction 
KB3 [18,19] produces a quasiparticle gap in 56 Ni about 1 MeV too large. Ap- 
proaching the doubly magic closure, the effects of this default become more 
visible. That is the reason why, in a recent study of the deformed excited band 
of 56 Ni [9], we used a preliminary modified version of KB3 in order to have 
the correct gap. This modified version of KB3 was also used in the study of 
the Ml strength functions of the N=28 isotones in ref. [20]. We shall exam- 
ine the interaction issues in section II. In sections III, IV and V we present 
the spectroscopic results for A=50, 51 and 52 respectively. In section VI we 
gather the beta decay results. In section VII we discuss the Coulomb Energy 
Differences (CED) between the yrast states of the mirror pair A=51. Finally, 
in section VIII we study the behaviour of the magnetic moments of the N=28 
isotones. We close with the conclusions. 

Throughout the paper / stands for f 7 / 2 (except, of course, when we speak of 
the pf shell) and r, generically, for any or all of the other subshells (pi/ 2 P3/2 fb/2)- 
Spaces of the type 

jn-no^no _|_ yn-no-l^no+1 _|_ . . . _|_ jn-no-t^no+t 

represent possible truncations: uq is different from zero if more than 8 neutrons 
(or protons) are present and when t = n — n we have the full space (pf) n for 
A = 40 + n. 
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The interaction KB3 is a (mostly) monopole modification of the original Kuo- 
Brown one [18]. The modifications are described in detail in [3]. 

In what follows, and unless specified otherwise, we use 

• harmonic oscillator wave functions with b = 1.01 A 1 / 6 fm; 

• bare electromagnetic factors in Ml transitions; effective charges of 1.5 e for 
protons and 0.5 e for neutrons in the electric quadrupole transitions and 
moments; 

• Gamow- Teller (GT) strength defined through 

B(GT) = (<£\ {„)', W = aW, (2) 
\9v / cS V ZJ i + 1 

where the matrix element is reduced with respect to the spin operator only 
(Racah convention [21]), ± refers to ft* 1 decay, t± = (t x ± 'iT y )/2, with 
t+p = n and {qa/ 9v)es is the effective axial to vector ratio for GT decays, 

( 9 -) =0.77(^) , (3) 

with (WtfvObarc = 1.2599(25) [22]; 

• for Fermi decays we have 

B(F) = (r) 2 , (r) = ^gffi (4) 

• half-lives, Ty 2 , are found through 

(f _l f e \ 1 1 6146 ±6 

Ua + I) 1/2 ~ Uv/fA)B(F) + B(GTY (5) 

We follow ref. [23] in the calculation of the Ja and f v integrals and ref. [24] 
for f e . The experimental energies are used. 

• The intrinsic quadrupole moments Q are extracted from the spectroscopic 
ones through 

go= 3^-j(J+l) gspcc(J) ' hlK ^ 1 ( 6 ) 
or from the BE2's through the rotational model prescription 

B (E2,J - J-2) = ^-e 2 \(JK20\J-2,K)\ 2 Ql for K ± 1 1. (7) 

The diagonalizations are performed in the m-scheme using a fast implementa- 
tion of the Lanczos algorithm through the code antoine [25] or in J-coupled 
scheme using the code NATHAN [26]. Some details may be found in ref. [27]. 
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The strength functions are obtained using Whitehead's prescription [28], ex- 
plained and illustrated in refs. [29-31]. 

All the experimental results for which no explicit credit is given come from 
the electronic version of Nuclear Data Sheets compiled by Burrows [32] . 



2 The Interactions 

Following ref. [33], we shall treat the effective interaction as a sum of a 
monopole and a multipole part. The monopole part is responsible for the 
energies of the closed shells (CS) and the closed shells plus or minus one par- 
ticle (CS±1). In our valence space, the starting point (or vacuum) is the 40 Ca 
core, and the single particle energies are provided by the levels of 41 Ca. The 
harmonic oscillator closure should be 80 Zr and its corresponding hole states 
in 79 Y. Nevertheless, and taking care of the effect of correlations, 48 Ca and 
56 Ni can be also taken as reference closed shells. In this we are lucky, because 
the information given by 80 Zr is rather useless. For N,Z>32 the influence of 
the lgg/2 orbit is very strong and the pf valence space is no longer valid. 
Around 80 Zr, the occupation of the orbits lgg/2 and 2d$/2 drives the nuclei 
into deformed shapes [34,35]. There are no experimental results available for 
the single hole states in 79 Y. Even if some excited levels were accessible, a 
particle plus rotor spectrum should be expected, coming from the coupling of 
the holes to the 80 Zr deformed core. With this information missing, we have 
to rely on indirect indications, as those coming from spherical Hartree Fock 
calculations (not very accurate) or from the new monopole formulae fitted to 
the real shell closures and their particle and hole partners [36]. Hence, our 
main guidance to adjust the monopole part of our interactions comes from 
the gaps around 48 Ca and 56 Ni and from the single particle states of 49 Ca and 
57 Ni. The multipole part of the interactions will be also analysed in terms of 
"coherent" multipoles, as proposed in ref. [33]. This part of the interaction is 
left unchanged. 

The quasiparticle gap for 48 Ca is defined as 

A = 2BE( 48 C&) - BE( 49 Ca) - 5£( 47 Ca) (8) 

and similarly for 56 Ni. 

In table 1 we compare these quantities for the most popular effective in- 
teractions used in the p/-shell; KB3 [18,19] and FPD6 [37]. Both interac- 
tions use very similar single particle energies; e7/2=0.0 MeV, €3/2=2.0 MeV, 
e V 2=4.0MeV and e 5/2 =6.5MeV for KB3 and e 7/2 =0.0MeV, e 3/2 =1.89MeV, 
e V 2=3.91 MeV and e 5 / 2 =6.49MeV for FPD6. The two body matrix elements 
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Table 1 

48 Ca and 56 Ni gaps and single particle energies, with the interactions KB3 and 
FPD6 (energies in MeV). 







A 






€5 — €3 

2 2 






e i —e3 

2 2 




A=48 


Exp. 


KB3 


FPD6 


Exp. 


KB3 


FPD6 


Exp. 


KB3 


FPD6 


t=0 


4.80 


5.25 


4.61 


3.59 


3.53 


2.66 


2.02 


1.70 


2.51 


Full 


4.80 


5.17 


4.69 


3.59 


3.80 


2.76 


2.02 


1.81 


2.37 


A=56 


Exp. 


KB3 


FPD6 


Exp. 


KB3 


FPD6 


Exp. 


KB3 


FPD6 


t=0 


6.39 


8.57 


7.41 


0.77 


0.38 


-0.48 


1.11 


1.15 


2.58 


t=3 


6.39 


7.73 


6.41 


0.77 


0.76 


0.07 


1.11 


1.14 


1.88 



defining FPD6 are scaled with the mass number A as (42/A) - 35 , while KB3 
does not incorporate any mass dependence. 



Notice that KB3 does definitely better than FPD6 for the single particle spec- 
tra of 49 Ca and 57 Ni, while for the 56 Ni gap, FPD6 is better. 

In addition the modification of KB3's gap, we want to to make it mass depen- 
dent, in order to be able to use it safely around and beyond 56 Ni, therefore we 
have to adjust the monopoles anew. As the gaps are subjet to variation when 
correlations are allowed, some trial and error fitting of the monopole changes 
is needed. The final modifications of KB3, defining KB3G for A=42 (we stick 
to the (42/A) 1 / 3 mass dependence) are the following: 



Vj= l (KB3G) = Vj= l (KB3) - 50keV, 
Vf=° (KB3G) = Vf= Q (KB3) - 100 keV, 
yjft (KB3G) = Vjff 2 (KB3) - 100 keV, 
V Ift (KB3G) = Vfff 2 (KB3) - 150 keV, 
VLp(KB3G) = V^(KB3) +400keV, 



where p denotes any of the orbits 2pi and 2p|. 

The T=0 and T=l modifications are different in order to recover simultane- 
ously the good gaps around 48 Ca and 56 Ni. The modification of the pp centroids 
bears no relationship with the gaps, it is aimed to give a single hole spectrum 
of 80 Zr in accord with the predictions of ref [36]. It is important to remark that, 
for the nuclei already studied with KB3, KB3G produces equivalent results. 
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Table 2 

Same as in table 1 with the modified interactions KB3G and FPD6*. 







A 






es — €3 






ei — e-i 




A=48 


Exp. 


KB3G 


FPD6* 


Exp. 


KB3G 


FPD6* 


Exp. 


KB3G 


FPD6* 


t=0 


4.80 


4.73 


4.61 


3.59 


3.20 


3.04 


2.02 


1.71 


2.13 


Full 


4.80 


4.69 


4.68 


3.59 


3.44 


3.10 


2.02 


1.82 


2.03 


A=56 


Exp. 


KB3G 


FPD6* 


Exp. 


KB3G 


FPD6* 


Exp. 


KB3G 


FPD6* 


t=0 


6.39 


7.12 


7.41 


0.77 


0.05 


0.24 


1.11 


1.23 


1.86 


t=3 


6.39 


6.40 


6.45 


0.77 


0.43 


0.68 


1.11 


1.19 


1.45 



In the case of FPD6, the gaps are nearly correct and the monopole defects 
amount to having the lf 5 / 2 orbit too low and the 2pi/ 2 orbit too high, both in 
49 Ca and 57 Ni. The modifications that repair that are: 



V Jft l (FPD6*) = Vjff; 1 (FPD6) + 50keV, 
Vf=^ (FPD6*) = Vf=^ (FPD6) - 50keV. 

In table 2 we present the resulting values around 48 Ca and 56 Ni. Notice that 
we have not aimed to an exact agreement. In particular, we know [39] that a 
t=3 calculations is not fully converged in 57 Ni, therefore we have let room for 
the H5/2 to move up and for the 2pi/ 2 to move down a little. We also show 
the evolution of the quasiparticle neutron gaps through the N=28 isotones 
for both modified interactions in table 3, The gaps produced by FPD6 do not 
differ appreciably from those of FPD6*. The agreement with the experimental 
results is the best that can be reasonably aimed to. 

Once the monopole part under control, we can analyse the multipole hamil- 
tonian using the method of ref . [33] . It amounts to diagonalize the monopole- 
free two body matrix elements in the particle-particle representation or in the 
particle- hole representation (i.e. after a Racah transformation). In the for- 
mer case, the physically relevant terms are the isovector and isoscalar pairing 
terms, while in the latter all the multipoles will show up with different coherent 
strengths given by the lowest eigenvalues of the respective matrices. In table 4, 
we have listed these numbers for the two interactions and for the most im- 
portant multipoles. The differences are small; basically consist in FPD6 being 
5-10% more intense than KB3 in all the channels. It is however satisfying the 
closeness of the values for two effective interactions derived by very different 
procedures. Remember that for its multipole part, KB3 is just equivalent to 
the original Kuo and Brown G-matrix [18] including the "bubble" correction, 
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Table 3 

Theoretical quasiparticle neutron gaps of the N=28 isotones (in MeV) compared 
with experiment. For Z>23 we list t=3 results, for the rest, full p/-shell results 
(energies in MeV). 







A 




C A / \ 

o = A{exp) — 


A(th) 


AT 1 

Nucleus 


FPD6* 


EXP. 


KB3G 


FPD6* 


KB3G 


48 Ca 


4.68 


4.81 


4.70 


0.13 


0.11 


49 Sc 


4.05 


4.07 


3.99 


0.02 


0.08 


50 rpi 
±1 


4.DD 


4.5 7 


A A 1 

4.41 


n no 
-0.09 


0.1b 


51y 


3.61 


3.7 4 


3.53 


0.13 


0.21 


52 Cr 


4.03 


4.10 


3.76 


0.07 


0.34 


53 Mn 


3.20 


3.12 


3.00 


-0.08 


0.12 


54 Fe 


4.30 


4.08 


4.02 


-0.22 


0.06 


55 Co 


4.00 


4.01 


4.16 


0.01 


-0.15 


56 Ni 


6.45 


6.39 


6.41 


-0.06 


-0.02 


Table 4 
Strengths (in 


MeV) of the coherent terms of the multipole Hamiltonian. 




Interaction 


particle-particle 




particle-hole 






JT=01 


JT=10 


Ar=20 


Ar=40 


Ar=ll 


KB3 


-4.75 


-4.46 


-2.79 


-1.39 


+2.46 


FPD6 


-5.06 


-5.08 


-3.11 


-1.67 


+3.17 


GOGNY 


-4.07 


-5.74 


-3.23 


-1.77 


+2.46 



while FPD6 was obtained via a potential fit to selected energy levels in the 
p/-shell. This supports of the conclusions in [33] about the universality of the 
multipole shell model Hamiltonian. 



We have purposely left to the end the line labeled GOGNY. It comes from the 
same analysis applied to the p/-shell two body matrix elements obtained using 
the density dependent interaction of Gogny [40]. The calculation was carried 
out using the single particle wave function obtained in a spherical Hartree- 
Fock calculation for 48 Cr in the uniform filling approximation in ref. [41]. 
In spite of the rather hybrid approach, the most important terms are again 
very similar to those arising from a G- matrix or from a shell model fit. In 
particular the agreement for the quadrupole and the spin-isospin terms is 
excellent. When it comes to pairing, this way of looking to the interaction 
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can help to overcome some language barriers between mean field and shell 
model practitioners. As it is evident from the table, the Gogny interaction has 
essentially the same amount of isovector and isoscalar pairing than the realistic 
interactions. Therefore, it contains the right proton-neutron pairing. Thus, if 
there is something to blame in the mean field calculations for N=Z nuclei, it 
should be rather the mean field approximations and not the interaction itself. 

In what follows we shall use mainly the interaction KB3G. We will compare 
in some cases with KB3 in order to evaluate the importance of the changes 
made in the monopole behaviour. The comparison with FPD6 will show to 
what extent the residual differences between good behaved interactions can 
have spectroscopic consequences. We have not attempted to compute all the 
states not even to draw those we have computed. Those not shown here can 
be obtained on request to the authors. 



3 The isobars A=50 

We have carried out the full p/-shell calculations for all the isobars. The 
results for 50 Cr and 50 Mn have been already published in refs. [16,4,17]. The 
experimental data for which no specific credit is given are taken from ref. [32]. 

3.1 Spectroscopy of 50 Ca 

The experimental data are compared with the calculation in fig. 1. All the ex- 
perimental states are plotted. We also plot all the calculated ones up to 5 MeV. 
Beyond, only the yrast states are drawn. The scarcity of spin assignment pre- 
cludes a more detailed analysis. The 2 + excitation energy is well reproduced 
as well as the gap in the spectrum between 1 and 3 MeV of excitation energy. 
KB3 gives equivalent results. No experimental information on transitions is 
available. 



3.2 Spectroscopy of 50 Sc 

Figure 2 reflects the experimental situation and our calculated spectrum. We 
have plotted all the experimental states even if only the lowest ones have 
assigned spin. All the calculated states up to the first 7 + state are shown. 
Above it, only the yrast band. The ground state multiplet, (I/7/2) 9 , 2p 3 / 2 is well 
reproduced. The first excited state must be a 2 + because the multiplet cannot 
contain two 3 + 's. The triplet at ~2 MeV, belonging to the configurations 
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; (8*,(6*)) 
■ (6*,<8*)) 



i (4* & 5*) 



1 4* 
2 + 

r 



Exp. 



KB3G 



Fig. 1. Energy levels of Ca 



9 
8 
7 
6 

>5 

Z A 

a 4 

3 
2 

1 h 





3* 
: 2*,3* 



= 2 + ,6 + ,3 + + + 
i ^ + „ + 4 ,6 
5 ,0 



Exp. 



KB3G 



Fig. 2. Energy levels 50 Sc. 

(I/7/2) 9 , (I/5/2, 2P1/2) 1 , is also well given. Beyond, the level density increases 
rapidly and no spin assignments are available. Notice that the experimental 
density of states is nicely reproduced up to 4 MeV. 
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L o* 

Exp. FPD6 KB3G 

Fig. 3. Energy levels of 50 Ti. 
3. 3 Spectroscopy of 50 Ti 

50 Ti is the stable member of the isobaric multiplet and the best known exper- 
imentally. It is semi- magic and therefore any defects in the neutron gap would 
be evident in the spectrum. In fig 3 the experimental level scheme is compared 
with the full pf-shell results using the interactions KB3G and FPD6. Up to 4 
MeV all the states are plotted. Between 4MeV and 5.6 MeV only those with 
unambiguous experimental spin assignment and beyond, only the yrast. The 
agreement is impressive for the two interactions, perhaps with a bonus for 
KB3G. Notice the excellent reproduction of the second and third + states as 
well as the perfect location of the high spin states. The second + corresponds 
to the excitation of two neutrons to the 2^3 orbit, while the third has a more 
complex structure. 



Table 5 

Transitions in 50 Ti. 



B(E2) (e 2 fm 4 ) 


Exp. 


Th. 


2+^0+ 


58(9) 


88 


4+^2+ 


60(1) 


86 


6+^4+ 


34(1) 


41 



The quadrupole and the magnetic moments of the yrast J = 2 + and J = 6 + 
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states are known [42]. Their values, ji exp (2 + ) = 2.89(15)/xat [43], Q e x P (2 + ) = 
+8(16) e fm 2 and n eX p(Q + ) — +9.3(10)/^ agree quite well with our predic- 
tions; /i(2 + ) = +2.5/ijv, Q{2 + ) = +Qefm 2 and /i(6 + ) = +8.3//^. A more 
detailed discussion on the magnetic moments is given in section 8. Similarly 
the E2 transitions of the yrast J = 2 + ,4 + y 6 + states are well reproduced by 
the calculation (see table 5). 



3.4 Energy levels of V 



The experimental information is also very rich for this nucleus. Given the 
high level density, we have represented in fig. 4 all the states up to 1.6 MeV 
only. Above, just the high spins and a couple of + 's and l + 's that may have 
experimental counterparts. The calculation gives a good reproduction of the 
ground state quintuplet and locates correctly the high spin states 9 + , 10 + and 
11 + . The only discrepancy affects to the bunch of 1 + , 2 + and 3 + states at 
around 1.5 MeV that are placed 500 keV too low. 
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> 



1 - 



■ ("*) 

• (101 

' (0*) 

• 0*,1* 

■ (9*) 

0*,1* 

■ 4 W ,5 W 

• (8)* 



(7)*,4* 
5* 



3*,2* 



Exp. 



KB3G 



Fig. 4. Energy levels of V. 

The magnetic and the quadrupole moment of the J = 6 + ground state have 
been measured [42]. Their values (/i exp (6 + ) = +3.3456889(14)/^, Q exp (Q + ) = 
+20.9(40) e fm 2 ) are very well reproduced by the calculations (/i(6 + ) = +3.18/zjv, 
Q(Q + ) = +19.6 e fm 2 ). The electromagnetic transitions of the yrast states are 
also known (table 6). Notice the excellent agreement found for the -dominant- 
Mi transitions. For the AJ=1, E2 transitions the huge error bars make the 
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comparison meaningless. On the contrary, for the only AJ=2 transition mea- 
sured, the agreement is very good. 

Table 6 

Transitions in 50 V. 



Exp. Th. 



B(M1) 


(/&) 




7+ -> 6+ 


1 2(2) 


1.0 


8+^7+ 


0.3(1) 


0.2 


11+ -► 10+ 


0.9(3) 


1.1 


B(E2) 


(e 2 /m 4 ) 


(e 2 / m 4 ) 


7+^6+ 


o 7c -+1313 
°' °-875 


108 


8+^6+ 


98(44) 


74 


8+^7+ 


219±g| 


10 


11+ -► 10+ 


iuy -109 


29 



4 The isobars A=51 

^. 1 Spectroscopy of 51 Ca 

The experimental data become rarer as we go far from stability. In 51 Ca only 
the ground state J = | is assigned in the experimental scheme. It is what 
is expected in any reasonable calculation. This nucleus was also studied in 
ref. [44], using the interactions KB3 and FPD6. As no comparison with the 
data is possible, we can examine the predictions of the different interactions. 
In fig. 5 we compare KB3 and KB3G. As the gap is very similar for the two 
interactions in this nucleus, the differences must be due to the changes in the 
pp interaction. This is clearly seen in the figure. The four states 3/2~, l/2~, 
5/2~ and 3/2~ are dominantly p 3 states and remain unchanged between KB3 
and KB3G. However, the states 7/2~ and 5/2~ are swapped with a relative 
change of about 2MeV. This is due to their different structures; the 5/2~ 
is /yP 2 /5, thus, with KB3G it gains 800 keV relative to the f^p 3 states, the 
7/2~ has a structure fjp 4 and it looses about 1.2 MeV. Similar arguments 
explain the lowering of the bunch of states at about 5 MeV excitation energy. 
Therefore a better experimental spectrum will be of much help in refining 
these monopolar changes. 
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8 

7 h 

6 
5 



(1/2,3/2,5/2) 



3 - 
2 - 
1 




(1/2,3/2,5/2) 



Exp. 



11/2" 
9/2" 



5/2" 
3/2" 



KB3 



KB3G 



Fig. 5. Energy levels of 5 Ca. 
4-2 Spectroscopy of 51 Sc 



As in 51 Ca, this nucleus was studied in ref. [44]. Some numerical errors present 
in this reference were discussed and corrected in refs. [45-47]. Figure 6 shows 
all the known experimental states. The calculated ones up to 4 MeV energy are 
also shown. Our results using KB3G are slightly better than those obtained 
with KB3 and substantially better than those produced by FPD6 [45,47], 
specially up to 2 MeV of excitation energy where the spacing of the calculated 
levels follows closely that of the experiment. 



4-3 Spectroscopy of 51 Ti 



In figure 7 we have plotted all the levels up to 3 MeV and a few more up 
to 4 MeV which have spin and parity assignments. No high spin states are 
known beyond this energy. The agreement is quite good, exception made of a 
couple of doublets that come out inverted (J = | and J —\ at 1.5 MeV 
and the second J = | with the first J = y at around 2.5 MeV). The 
experimental spin assignment for the state around 3.7 MeV is either J = 
— or J = y . Our calculation predicts a J = y close to this energy, 

while the first J = y state is predicted 2 MeV higher. If we compute the 
51 Ti spectrum with FPD6 there are not significant differences, because the 
dominant configurations are (fip\)u (fz)n, therefore the /s does not play an 
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1 



Exp. 



■ (1/2,3/2,5/2) 
- (3/2",5/2") 

; (3/2-.5/21 

■ (3/2") 



: (1/21 



. (3/2,5/2) 

■ (3/2,5/2,7/2*) 
' (11/2)- 

■ (3/2)- 



(7/2)" 



KB3G 



: 5 /2- 9,2 ~ 
■ 11/2,3/2- 



11/2 
3/2- 



Fig. 6. Energy levels of 51 Sc. 



2 



(13/2,17/2) 



(7<2.S«2: 
^5/2,7/2)- 

(15/2-) 

(7/2,9/2)- 

7/2- 

(11/2") 
3/2" 



(5/2)- 



9/2" 
11/2- 



Exp. 



KB3G 



Fig. 7. Energy levels of Ti. 

important role. We can interpret most of the 51 Ti level scheme as the result 
of the coupling 50 Ti®p\ . 



Some electromagnetic transitions are experimentaly known that are well re- 
produced by the calculation (see table 7). Again in this case, the E2 transitions 
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with AJ=1 are poorly determined experimentally, due to the dominance of 
the Ml transitions. The calculated AJ=2 B(E2)'s agree extremely well with 
the measured values. 

Table 7 

Transitions in 51 Ti. 





Exp. 


Th. 


B(M1) 






5- 3- 
2 2 


0.06(2) 


0.07 


B(E2) 


(e 2 /m 4 ) 


(e 2 /m 4 ) 


5- 3- 
2 2 


348(112) 


88 


7- 3~ 
2 2 


225(202) 


88 


11- 7- 
2 2 


95(17) 


95 


15-11- 
2 2 


62(24) 


53 



^.^ Spectroscopy of V 

51 V is the stable isotope in the isobar chain A = 51 and corresponds to the 
N=28 neutron shell closure. The experimental information is very rich, ex- 
tending up to 13MeV. Figure 8 shows the experimental yrast band with all 
the known high spin assignments compared with the calculated one includ- 
ing three more predicted spins. The agreement is extremely good. For a more 
complete analysis below 3 MeV energy, where more experimental information 
is available, we have plotted in figure 9 all the levels up to that energy and 
a few more above that can be put in correspondence with the experiment. 
We have calculated with KB3G and FPD6, both reproduce nicely the level 
scheme but definitely the quality of the agreement with KB3G is better. The 
magnetic moments of the ground state J — \ and first excited state J = | 
as well as the quadrupole moment of the ground state are known [42]. Their 
values are: 

Vex P ( 7 f) = +5.14870573(18) /i N Q exp ( 7 f) = -5.2(10) e fm 2 
Hex P (l~) = +3.86(33) fx N 

while the calculation gives: 

//(|") = +4.99 QCf) = -Q-5efm 2 
Mf ~) = +3.36 pL N 
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Fig. 8. Yrast band of 51 V. 
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Fig. 9. Energy levels of 5 V. 

in good agreement with the data. Furthermore, some electromagnetic transi- 
tions, mainly between yrast states are known experimentally [42] . The calcu- 
lation does also quite well, as it is shown in table 8. 
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Table 8 

Transitions in 51 V. 



Exp. Th. 


B(M1) 






9- 7- 
2 2 


0.0006(2) 


0.2-10~ 6 


1 ^ — 1 1 — 

2 2 


<0.0077 


0.014 


B{E2) 


(e 2 /m 4 ) 


(e 2 /m 4 ) 


9- . 7- 
2 2 


35(6) 


32 


11- , 7- 
2 2 


95(8) 


103 


13- 11- 
2 2 


<34.8 


0.01 


15- 11- 
2 2 


66(6) 


78 



4-5 Spectroscopy of 51 Cr 



We present in fig. 10 the results for the yrast band of 51 Cr in a t = 5 trunca- 
tion and in the full space using KB3G. We have also included the low- lying 
J = \ , § and I states, and a second state of each spin from J = y 
on. The figure shows that the full calculation do not modify substantially the 
results obtained at t — 5. The agreement with the experimental data is very 
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Fig. 10. Yrast band of 51 Cr; experiment, t = 5 and full calculation. 
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satisfactory. Notice that our prediction for the spins of the triplet of states 
around 6 MeV is shifted down by one unit relative to the preliminary exper- 
imental assignment, that should be revised in the light of our results. The 
detailed comparison with all the levels up to 2.8 MeV, and with the high spin 
states above this energy up to J = y is carried out with a t = 5 calculation. 
The results for both KB3 and KB3G are plotted in figure 11. 

The agreement is very good in spite of some local inversions. We can also 
notice that KB3G improves systematically the KB3 results. The sequence 
of states J = | , | , J = y belong to the configuration (I/7) 11 , while 

the doublets J = \~ ,\~ and J = §~, \~ located at 800 keV and 1.5 MeV, 

correspond to a neutron jump to the orbit 2j93, a sort of coupling 50 Ci®p\. 

2 2 

As a consequence of the reduction of the gap, KB3G put these doublets at 
lower energy, while leaving the yrast ones unchanged. This provides a good 
illustration of a direct monopole effect. A similar though more pronounced 
effect is seen in the lowering of the doublet J = § , § experimentally at 
2 MeV predicted by KB3 500 keV too high, again a gap effect solved by 
KB3G. Besides, the high spin states were also too high in KB3 and come now 
very close to their experimental position. 
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Fig. 11. Energy levels of Cr, t = 5. 



The magnetic moments of the ground state and the first excited state are 
known, (/i exp (| ) = —0.934(5) and /i exp (| ) = —0.86(12) jj, N ). Our predic- 
tion for the ground state is quite accurate /i(| ) = —0.886^ while for the 
J = § state our result /i(| ) = —0.40 fi N is off by a factor two. 
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In table 9 we compare the experimental and calculated transition probabilities. 
The B(Ml)'s and the AJ=2 B(E2)'s are well reproduced . However this is not 
the case for the AJ=1 B(E2)'s, in particular the discrepancies are huge for 
the transitions -y — > | and ^ — > ^ ■ We have carried out a consistency 
test using the experimental M1/E2 branchings, 5, and we have found that 
the experimental 5 values are inconsistent with the experimental B(E2)'s but 
consistent with the calculated values. This is not surprising, because of the 
complete dominance of the Ml transition, that may cause large errors in the 
extraction of the B(E2) value. In order to discard any other origin of the 
discrepancy we have repeated the calculation with FPD6 and the situation is 
the same or even worse. 
Table 9 

Transitions in 51 Cr. 





Exp. 


Th. 


B(M1) 


H) 


(fit) 


9- , 7~ 
2 2 


0.31683(3043) 


0.439 


11- 9~ 
2 2 


1.253(537) 


1.295 


13-11- 
2 2 


0.8950(1969) 


1.356 


B(E2) 


(e 2 /m 4 ) 


(e 2 fm 4 ) 


9- _, 7- 
2 2 


124(56) 


213 


11- 7- 
2 2 


67(34) 


72 


11- 9- 
2 2 


8(3) 


180 


13- 11- 
2 2 


6(1) 


151 


15-11- 
2 2 


44(1) 


53 



4-6 Spectroscopy of 51 Mn 

As for 51 Cr, we present in fig. 12 the comparison between the experimental 
level scheme and our calculations in the full space and with a t = 5 truncation 
for the yrast band and the | , | doublet around 2 MeV. Again no relevant 
difference is appreciated between both calculated schemes and the agreement 
with the experiment is very good. Therefore, a t = 5 truncation is used to get 
a more detailed spectroscopy for all the states below 3 MeV and for the high 
spins up to J = f (fig. 13). New data are available from a Gammasphere 
experiment aimed to the study of Coulomb energy differences in mirror nu- 
clei [48]. The agreement is very good in both regions, although the high spin 
states are slightly shifted up. 
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Fig. 12. Yrast band of Mn; experiment, t = 5 and full calculation. 
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Fig. 13. Energy levels of 5i Mn, t = 5. 

Good agreement is also obtained for the electromagnetic moments of the 
ground state. The experimental values: 



f^exp\ 



Q 



exp\2 



= 3.5683(13) n N 
") = 42(7) efm 2 
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are very well reproduced by the KB3G calculation: 

= 3.40 n N 
Q(jj-) =35e/m 2 

In addition, there is experimental information on electromagnetic transitions, 
that we compare with the calculation in table 10. The accord is excellent for 
the Ml transitions. For the E2's the situation is similar to the other isobars; for 
the AJ=2 transitions the agreement is quite good, while some discrepancies 
are present in the AJ=1 cases, that can be very large, for instance in the 
transition from J = y to J = | . As in the 51 Cr case, we have found 
an inconsistency between the experimental values for the B(M1) and B(E2) 
and the experimental mixing parameter 5. On the contrary the experimental 
value is fully compatible with the computed transition probabilities. Notice the 
abrupt change (more than two orders of magnitude decrease) in the transition 
probabilities either Ml or E2 of the J — y state, that are spectacularly 
reproduced by the calculation. A very intuitive physical explanation of this 
isomerism will be given later, when discussing the Coulomb energy differences 
between the yrast states of the mirror pair 51 Fe- 51 Mn. 

Table 10 

Transitions in 51 Mn. 





Exp. 


Th. 


B(M1) 




04) 


7~ 5- 
2 2 


0.207(34) 


0.177 


9- . 7~ 
2 2 


0.16(5) 


0.114 


11- 9~ 
2 2 


0.6623(2148) 


0.423 


17- 15- 
2 2 


0.00012(4) 


0.00003 


19- 17- 
2 2 


>0.5728 


0.801 


B{E2) 


(e 2 /m 4 ) 


(e 2 /m 4 ) 


7- 5- 
2 2 


528(146) 


305 


9- 5- 
2 2 


169(67) 


84 


9" 7~ 
2 2 


303(112) 


204 


11- 7~ 
2 2 


236(67) 


154 


11- 9~ 
2 2 


4.16(135) 


190 


17- 13- 
2 2 


1.236(337) 


2.215 
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5 The isobars A=52 



As in the A = 51 isobaric multiplet we have performed full Oku calculations 
for all the isotopes in 52 Cr and 52 Mn the full calculation is limited to the yrast 
levels. More detailed spectroscopy is carried out with a t—5 truncation. The 
results of the full p /-shell calculation for 52 Fe have been already published 
in [5]. The experimental values are taken from ref. [32]. 



5.1 Spectroscopy of 52 Ca 



The experimental information is scarce, as can be seen in figure 14. Both KB3 
and KB3G give the same excitation energy for the doublet J = l + ,2 + at 
2.5 MeV, and reproduce correctly the energy of the first excited state J = 
2 + . The states whose leading configuration has four p-particles behave like 
the ground state under the KB3G changes, therefore we expect that their 
excitation energies do not change between KB3 and KB3G. This is the case 
for the above mentioned doublet, whose leading configuration is /f p\p\ and 

2 2 2 

for the second + whose configuration is /f pfpl, therefore, both interactions 

2 2 2 

locate them at the same place. 
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Fig. 14. Energy levels of 52 Ca. 

The multiplet J = l + ,2 + ,3+,4+ predicted at 5 MeV by KB3 has a configu- 
ration /f p I fl , consequently it is less affected by the monopole corrections 
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(AVpp) made in KB3G than the ground state, and moves down 1 MeV in ex- 
citation energy. The states J = 5 + , J = 6^, J = 7 + , J = 8 + and J = 9 + are 
also lowered by KB3G, due to the occupation of the \j\ orbit. 



5.2 Spectroscopy of 52 Sc 



The dominant configuration in the ground state of 52 Sc is (fz)ir(fzp\)v It 

2 2 2 

produces the multiplet of states with J going from 2 + to 5 + , below 1 MeV. 
Experimentally little is known, the ground state 3 + , one state without spin 
assigned at about 0.6 MeV and four 1 + states seen in the beta decay of 52 Ca [49] 
(see fig. 15). In the figure we include the calculated ground state multiplet and 
the high spins up to J = 8 + , as well as all the J = 1 + states present below 
4.5 MeV. According to the KB3G results, the state at 0.6 MeV should be the 
2 + member of the ground state multiplet. In ref. [50] it is argued that this state 
could be a 1 + , based on the FPD6 results. However this has to be attributed 
to the fact that FPD6 puts the l/s orbit too low at N=28. If this state were 
a 1 + , it would have been seen in the decay of 52 Ca, and this is not the case. 
The calculation produces six 1 + states in the region where experimentally only 
four have been found. We will discuss the decay pattern of 52 Ca in section 6.3. 
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Fig. 15. Energy levels of 52 Sc. 
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5. 3 Spectroscopy of 52 Ti 



The experimental level scheme is compared in figure 16 with the full pf-shell 
results using the interactions KB3G and FPD6. The location of the yrast 
states from J = + to J = 6 + corresponding to the dominant configuration 
(fi Pi)v(fi)n i s correct for both interactions. The lowest states of 52 Ti can 

2 2 2 

be described essentially in terms of the coupling 50 Ti(g>pf . Notice the right 

2 

position of the triplet J = 2 + , 4 + , 2 + at ~ 2.4 MeV. However, the doublet 
J = 4+ ,2+ at ~ 3.5 MeV is clearly better placed by KB3G than by FPD6. 
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Fig. 16. Energy levels of 52 Ti. 



These states are dominantly /f/f and they are shifted down by FPD6. This 

2 2 

is even more evident for the second J = + state that is reasonably positioned 
by KB3G but strongly shifted down by FPD6. The richness of the calculated 
level scheme calls for an improvement of the experimental spectrum. 



There are six experimental transitions known [42], though with large uncer- 
tainties. The calculated values are in accord with the data, as can be seen in 
table 11. 
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Table 11 

Transitions in 52 Ti. 



Exp. Th. 

B(Ml) (^) 

2+ - 2+ 0.551°;^ 0.85 

2^ 2^ >0.16 0.51 

B(£?2) (e 2 /m 4 ) (e 2 /m 4 ) 

2^ -► 0+ 138l^ 4 85 
2$ -► 0+ 31±f| 16 

2^ -► 2^ >127 66 

6+ -► 4+ 123(22) 80 



5.4 Spectroscopy of 52 V 



In figure 17 we have plotted all the states with parity and spin assignments up 
to 2 MeV. Above this energy only the levels with spins equal or greater than 
6 + are shown. The situation is similar to 52 Sc. The dominant configuration 

(fi)n(fip\)v in the ground state produces the multiplet of states with J 

2 22 

going from 1 + to 5 + . 
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Fig. 17. Energy levels of 52 V. 
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The calculated ground state multiplet is more dilated that its experimental 
counterpart. Although the more relevant features are well accounted for we 
do not predict the correct ground state spin. The J = 3 + ,4 + doublet around 
850 keV comes out inverted. These states can be described as the result of the 
coupling 51 V®£>3. The high spin states J = 7 + and J = 9 + are shifted up, 

2 

making the agreement with experiment worse than in the other nuclei studied 
in this work. 

Only two electromagnetic transitions are experimentally known [42]: the E2 
connecting the states J = 9+ and J = 7+ (B(E2) = 81(6) e 2 fm 4 ) and 
the Ml conecting the first excited J = 2 + state and the J = 3 + ground 
state (B(M1) = 1.7(4) 

A* at)- F° r the first one the calculation produces a value 
B(E2) = 87 e 2 fm A and for the second £(M1,2+ -> 3+) = 1.7/2%, a rather 
impresive agreement. 

5.5 Spectroscopy of 52 Cr 

Figure 18 shows the yrast states of 52 Cr, up to the band termination, calculated 
in the full p/-shell and in a t = 5 truncation, compared with the experimental 
values. It is quite evident that the states obtained with a t = 5 truncation 
are nearly converged. The accord of the KB3G results with the experiment is 
excellent. 
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Fig. 18. Yrast band of 52 Cr; experiment, t = 5 and full calculation. 
More detailed spectroscopic results are shown in figures 19 and 20 in a t = 5 
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truncated calculation. In fig. 19 all the states experimentally known up to 
~ 4 MeV [32] have been included as well as the corresponding ones obtained 
with KB3G and KB3. The levels above this energy are plotted in figure 20, 
where the experimental data come from a recent experiment [51]. The agree- 
ment for the non yrast states is also very satisfactory. Notice the good po- 
sitioning of the doublets J = + ,4 + at 2.7 MeV and J = 4+,3 + at 3.5 MeV 
and the one to one correspondence between the three J = 2 + and the three 
J = 5 + states measured in the excitation energy range 3 MeV~ 4 MeV and 
the calculated ones. Beyond 4 MeV the yrare states and the odd yrast are also 
fairly well accounted for. It is clear from fig. 20 that KB3G produces a more 
satisfactory level scheme than KB3, that shifts up all the levels with spins 
greater than J = 8 + , the maximum value attainable within the (I/7/2) 12 con- 
figuration. This is again a manifestation of the too large N=28 gap produced 
by this interaction. 
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Fig. 19. Energy levels of 52 Cr up to 4 MeV, t = 5. 

The experimental magnetic and quadrupole moments for the first excited state 
J = 2+ are: 



Qexpi^ ) 



+3.00(50) fi N 
-8.2(16) e/m 2 



and the calculated ones: 



/i(2+) = +2.50 //jv 
Q(2+) = -9.4 e/m 2 
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Fig. 20. Energy levels of 52 Cr above 4MeV, t = 5. 

In table 12 the known experimental values for the electromagnetic transitions 
between the states of the yrast sequence are given. The correspondence with 
the predicted ones is good with the exception of the E2 transition J = 4 + to 
J = 2 + , whose measured value is far outside the range expected in this mass 
region. Given its very large error bar, we would rather disregard this measure. 



Table 12 

Transitions in 52 Cr. 





Exp. 


Th. 


B(M1) 


(fit) 


04) 


9+^8+ 


0.05728(3759) 


0.040 


B(E2) 


(e 2 /to 4 ) 


(e 2 /to 4 ) 


2+^0+ 


131(6) 


132 


3+^2+ 




5 


4+^2+ 


761(265) 


107 


6+^4+ 


59(2) 


68 


8+^6+ 


75(24) 


84 


9+^8+ 


0.5(20) 


0.6 
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5.6 Spectroscopy of 52 Mn 



The full space calculation for the yrast states of the odd-odd nucleus 52 Mn 
produces only minor differences with the t — 5 results (fig. 21). Notice the full 
correspondence between the states of the multiplet below 1 MeV. The exper- 
imental data for the spins beyond the band termination (11 + to 16 + ) come 
from a recent experiment [52]. The agreement between experiment and theory 
is spectacular. There is only three experimental states with spin assignment 
not represented in the figure: a second J = (5 + ) at 1.42 MeV, a third J = (5 + ) 
at 1.68 MeV and a second J = (6 + ) at 1.96 MeV. The t — 5 calculation places 
the second J = 5+ at 1.47 MeV, the third J = 5+ at 2.05 MeV and the second 
J = 6+ at 1.91 MeV. 
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Fig. 21. Yrast band of 5 Mn; experiment, t = 5 and full calculation. 



The experimental magnetic and quadrupole moments for the ground state are 
known: 



Ve Xp (6 + ) = +3.0622(12) fi N 
Q exp (6 + ) = +50(7) efm 2 

the calculated values reproduce nicely these values: 

^(6+) = +2.9518 fx N 
Q(6+) = +50e/m 2 
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Some electromagnetic transitions along the yrast sequence have been mea- 
sured. The calculations are in fair agreement with them (see table 13). 

Table 13 

Transitions in 52 Mn. 





Exp. 


Th. 


B(Ml) 




04) 


7+^6+ 


0.5012(2506) 


0.667 


8+ -► 7+ 


>0.015931 


0.405 


9+ -> 8+ 

U 7 O 


1 074+H^ 

±.Ui<i_ Q 537 


759 


B(E2) 


(e 2 fm 4 ) 


(e 2 fm 4 ) 


7+^6+ 


Q9 +484 
yz_ gl 


126 


8+^6+ 


>1.15 


33 


8+^7+ 


>4.15 


126 


9+^7+ 


104+joo 


66 


11+ 9+ 


54(6) 


53 



6 Half-lives and other /3-decay properties 

Once the level schemes have been analysed, we study the f3~ decays of these 
nuclei, to complete the description of this mass region. We present the results 
with the same ordering used for the level schemes. Because of the increasing 
sizes of the calculations we have limited this study to the isotopes of Calcium, 
Scandium, Titanium and Vanadium. The calculated values for KB3 and KB3G 
gathered in the tables correspond to full Ohw calculation unless otherwise indi- 
cated. The Qp- values are also included. The errors attached to the calculated 
values proceed from the errors in the experimental Qp- values. We compute 
the half-lives by making the convolution of the strength function produced 
by the Lanczos method with the Fermi function, increasing the number of 
iterations until convergence is achieved. 

6.1 j3~ decays in the isobar chain A = 50 

The experimental and calculated half-lives in the isobar chain A = 50 are 
given in table 14. The agreement is quite satisfactory. 
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Table 14 



Half-lives of A 


= 50 isobars. Q^- 


values from [32]. 










Ti 

2 






a r 


KB3 


KB3G 


Exp. 


Q™ p (MeV) 


50 C a 0+ 


11 l +0 - 3 s 

11.1_ 2 6 


11.210.3 s 
12.3 ± 0.3 s 


13.9(6) s 


4.966(17) 


50 Sc 5+ 


130±3s 


1401? s 
120+1 s 


102.5(5) s 


6.888(16) 


50 Ti 




Stable 




50 y 




4 th forbidden C.E. 





The Gamow- Teller strength function for 50 Ca is shown in fig. 22. The amount 
of strength below the Qp- value is very small and it is concentrated at ~ 
2MeV. Notice that there is a large amount of strength close above the Qg- 
threshold. Thus, small variations in the excitation energies of the daughter 
states could influence (although not drastically) the final value of the half-life. 
In 50 Sc (see figure 23) the situation is similar. Only a tiny fraction of the total 
strength lies below the Qp- value. 



6 - 
5 - 




5 10 15 20 25 30 

E(MeV) 

Fig. 22. Gamow- Teller strength after 50 iterations for the 50 Ca 0~ decay. The 
< ot > 2 values are summed up in 1 MeV bins. The dashed line indicates the exper- 
imental Qp- value. 

We have computed the intensity (in percentage) of the ground state decay to 
the levels within the Qp- window of the daughter nucleus. 50 Ca decays 99.7% 
to the excited state J = 1 + at 1.78 MeV in 50 Sc, what agrees perfectly with the 
experimental situation (99.0(13)% decay to the J = 1 + state at 1.85 MeV). 
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6 



5 - 



4 - 



2 




Fig. 23. Gamow- Teller strength after 50 iterations for the Sc (5~ decay. The 
< err > 2 values are summed up in 1 MeV bins. The dashed line indicates the exper- 
imental Qp- value. 

Analogously, figure 24 shows the experimental and calculated intensities for 
the f3~ decay of 50 Sc. Both experimental peaks are reproduced by the calcula- 
tion although they are slightly shifted. Since the 50 Ca decay feeds an unique 
state and 50 Sc decays to just two, it is possible to use the experimental en- 
ergies of the states to calculate the half-lives and to gauge the effect of the 
phase space. Proceeding in such a way, the new calculated half-lives for 50 Ca 
and 50 Sc are, respectively, Ti = 12.310.3 s and Ti = 1201? s which improve 
the agreement with the experimental half-lives. These results are also shown 
in table 14. 
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Fig. 24. Percentages of the j3 decay of 50 Sc. Experimental data from [32]. 
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6.2 (3 decays in the isobar chain A = 51 



Table 15 gathers the experimental and calculated half-lives for the isobar chain 
A = 51. The calculated values are quite close to the experimental ones. 

Table 15 

Half-lives of A = 51 isobars. Qp- values from [32]. 









Ti 

2 






A 


r 


KB3 


KB3G 


Exp. 


QjI(MeV) 


51 Ca 


2 




7.618;?. 


10.0(8) s 


7.332(93) 


51 Sc 


7- 

2 


14.110.3 s 


10.410.2 s 
12.3 ± 0.3 s 


12.4(1) s 


6.508(20) 


51 Ti 


3- 
2 


6.9410.02 mm 


8.04l0.02min 


5.76(1) rain 


2.4706(15) 



6.97 ±0.02 min 
V Stable 



Although within the Qp- energy window there is a large amount of Gamow- 
Teller strength (see figures 25, 26 and 27), the situation looks like in the A = 50 
chain and much strenght is located around the Qp- values. 




Fig. 25. Gamow- Teller strength after 50 iterations for the Ca (3~ decay. The 
< ot > 2 values are summed up in 1 MeV bins. The dashed line indicates the exper- 
imental Qp- value. 

In figures 28 and 29 we compare the experimental percentages for the decays 
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5 10 15 20 25 30 35 40 
E(MeV) 



Fig. 26. Gamow- Teller strength after 50 iterations for the Sc (5~ decay. The 
< err > 2 values are summed up in 1 MeV bins. The dashed line indicates the exper- 
imental Qp- value. 




Fig. 27. Gamow- Teller strength after 50 iterations for the Ti (3~ decay. The 
< ot > 2 values are summed up in 1 MeV bins. The dashed line indicates the exper- 
imental Qp- value. 

with the calculation for both 51 Ca and 51 Sc. In the latter, the agreement is 
very good. 51 Ca demands some comments. To compare with the experimental 
data we must keep in mind that only states up to 3.8 MeV in 51 Sc have been 
observed. In the first MeV there is a single experimental state fed by the decay. 
This fact is well reproduced by the calculation. Up to 3.5 MeV there are seven 
states experimentally observed in the decay, the same number is given by the 
calculation. However, the distribution of the intensity is not the same. In the 
calculation all the intensity is concentrated between 2.5 MeV and 3.0 MeV, 
while experimentally there are two states strongly fed at lower energy. 
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Fig. 28. Percentages of the f3 decay of 51 Ca. Experimental data from [32]. 
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Fig. 29. Percentages of the j3 decay of 51 Sc. Experimental data from [32]. 

In the case of 51 Ti the decay only feeds two states in the daughter nucleus, 
a J = |" at 0.928 MeV (8.1(4)%) and a J = §~ at 0.320 MeV (91.9(4)%). 

The calculation agrees with it since it predicts the decay to a J = | state at 

1.155 MeV (6.5%) and to another J = §" state at 0.457 MeV (93.5%). 

In view of the comparison of the intensities for the decays of 51 Ca, 51 Sc and 
51 Ti, it seems plausible to make a new analysis in the same terms we did in 
the previous section, i.e. to take the experimental energies for the four states 
strongly fed in 51 Sc and the corresponding two in 51 Ti. The new half-lives we 
obtain are Ti = 12.310.3 s for 51 Sc and Ti = 6.97+0.02 mm for 51 Ti, in much 

2 2 

better agreement with the experimental values shown in table 15. For 51 Ca 
this analysis is not so easy since the number of states fed is larger and it is 
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not simple to make a one to one correspondence between experimental and 
theoretical states. 

6.3 f3~ decays in the isobar chain A = 52 

In table 16 the experimental data and the calculated results for the half-lives 
of the nuclei in the A = 52 isobar chain are shown. For 52 V we have only made 
a t = 4 calculation, because, to be consistent, the calculation in the daugther 
has to be t = 5. Going beyond that would have demanded an enormous 
computational effort. The theoretical half-lives compare reasonably well with 
the experimental numbers, although the discrepancy in 52 Ca is severe. We 
shall now analyse the intensity distributions and try to recompute the half- 
lives using as much as possible the experimental excitation energies. 

Table 16 

Half-lives of A = 52 isobars. Qp- values from [32]. 









Ti 

2 






A 


r 


KB3 


KB3G 


Exp. 


QjI(MeV) 


52 Ca 


0+ 




1.38l°;Is 


4.6(3) s 


7.900(500) 


52 Sc 


3+ 


8t 2 lS 


6 2 +1 -° s 


8.2(2) s 


9.010(160) 


52 Ti 


0+ 


2.6910.04 min 


3.3810.07 mm 
2.30 ± 0.04 mm 


1.7(1) min 


1.973(8) 


52y 


3+ 


5.544±o$4 min 


5.7910.01 min 
4.87 ±0.01 min 


3.743(5) min 


3.9756(12) 



Figure 30 shows the percentages for the Ca (3~ decay and figure 31 the same 
for 52 Sc. The 52 Ca case demands a special attention. When we analyzed the 
52 Sc spectrum in section 5.2 we mentioned that there was one state without 
any experimental spin assignment at 675 keV (fig. 15). In [50] the calculation 
was made with the FPD6 interaction, producing a first J = 1 + at 701 keV. It 
was suggested that it could correspond to the 675 keV state. Our interpretation 
of this state as the 2 + member of the ground state multiplet is experimentally 
supported by the fact that if this state were a J = 1 + , it would have been 
observed in the f3~ decay of 52 Ca. This is not the case, as can been seen in 
figure 30. On the contrary, most of the decay (86.8%) goes to the first exper- 
imental state with J = 1 + assignment, corresponding to the one we obtain in 
the calculation. Thus, it is discarded that this unassigned experimental state 
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Fig. 30. Percentages of the f3 decay of 52 Ca. Experimental data from [32]. 
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Fig. 31. Percentages of the f3 decay of 52 Sc. Experimental data from [32]. 

could be a J = 1 + . The first J = 1 + corresponds to a configuration with a 
neutron in the fs orbit. As we said in the introduction, FPD6 locates the lfs 

J 2 ' J 2 

orbit too low, producing a very low J = 1 + state. The other three experi- 
mental J = 1 + states are also fed through the decay of 52 Ca although with 
smaller intensity. The percentage of the second at 2.7 MeV is 11.2% while for 
the third at 3.5 MeV it is 0.6% and 1.4% for the last one at 4.3 MeV. These 
data help us to make a one to one correspondence between the states above 
2 MeV in figure 15. Our third state J = 1 + at 2.6 MeV corresponds to the 
second experimental one, our fourth J = 1 + at 3.4 MeV to the third and the 
fourth experimental J = 1 + to our sixth at 4.4 MeV. 

The situation for 52 Sc in figure 31 is very satisfactory. The decay to the J = 2 + 
excited state at lMeV and to the triplet around 2.3 MeV is well reproduced, 
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although in the calculation the triplet is more expanded than the experimen- 
tal one. There is one state at 3.1 MeV without assigned spin, that seems to 
correspond to the 3 + state near 3 MeV (the calculation predicts a 0.7% branch 
to it). The next branch goes to the second 4 + at 3.3 MeV. Beyond, it becomes 
difficult to establish a detailed correspondence. 

The decay of 52 Ti feeds (100%) a single state in the daughter nucleus 52 V, 
the J = 1 + at 141 keV. This fact is reproduced properly by the calculation 
which predicts a 100% feeding of the J = 1 + state at 300 keV in the daughter 
nucleus for the (3~ decay of 52 Ti. On the other hand, 52 V decays 99.22(5)% 
to the J = 2 + excited state of 52 Cr located at 1.4 MeV, what agrees with the 
calculation (a 98.92% feeds the J = 2+ state at 1.5 MeV). 

Finally, we recalculate the half-lives using the experimental energies of the 
single state fed in the decay of 52 Ti and 52 V, and of the two states fed in the 
decay of 52 Ca. For 52 Sc the analysis in these terms is not so easy. The new 
calculated half-lives are Ti = 2.3010.04 mm for 52 Ti and Ti = 4.8710.01 min 

2 2 

for 52 V. Both results improve the previous ones shown in table 16. For 52 Ca we 
get a half-life Ti = 1. 381^4 s, reducing substantially the initial discrepancy. 



7 Coulomb energy differences (CED) in the mirror pair 51 Mn- 51 Fe 

It is well known that if the nucleon-nucleon interaction were charge symmetric 
the mirror nuclei would have the same level scheme. Thus, the small differences 
(normally a few tens of keV) in their excitation energies, are due to the isospin 
symmetry breaking Coulomb interaction. The analysis of this effect has been 
recently pushed up in mass, with the study of the A=47 and A=49 mirror pairs 
[53,54]. It was found that the CED's are extremely sensitive to the structure 
of the nuclear wave functions. The next step has been to measure the 51 Mn- 
51 Fe mirror pair, the heaviest one in which high spin states up to the band 
termination have been observed [48]. 

In figure 32 the calculated spectra for both nuclei, including the Coulomb 
interaction are compared with the experimental data. The Coulomb matrix 
elements are those denoted "A42" in ref [54]. At this scale, very slight dif- 
ferences (hardly any) are appreciable between them. However, looking more 
carefully at the CED's an abrupt change is observed at J = y (fig. 33). The 
effect is also present in the calculation, that shows exactly the same trends as 
the experiment, although with enhanced values. The large increase in the CED 
can be interpreted as due to the alignment of one proton pair in 51 Fe which 
does not occur in 51 Mn. As a consequence, the Coulomb energy is sharply 
reduced in 51 Fe but not in 51 Mn, therefore the CDE increases dramatically (in 
absolute value). The breaking of a proton pair occurs in 51 Fe because it has 
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Fig. 32. Excitation energies of yrast bands of the mirror pair Mn- Fe. KB3G in 
the full p/-shell. Coulomb included. 

an even number of protons while the odd proton in 51 Mn blocks the break- 
ing process. In 51 Mn, a neutron pair is broken, but this has no effect in the 
Coulomb energy. Beyond J = y~, the protons start aligning also in 51 Mn and 
therefore the CED's approach zero at the band termination. 




Fig. 33. Experimental and calculated CED defined as E x ( 51 Fe) - £ x ( 51 Mn) 



In order to make more visible the proton pair breaking, we have calculated 

-|0 



the expectation value of the operator H a [ in 



\6,1 / 



\6,1 



for the 



states of each nucleus. This operator acts as a sort of counter for the pairs of 
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l/i protons coupled to the maximum spin value J = 6. In figure 34 we have 
plotted the difference between its expectation values for the states of 51 Mn 
and 51 Fe. It is manifest that a proton pair fully aligns at J = y in 51 Fe and 
not at all in 51 Mn. 




-1.0 1 — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — 1 

3 5 7 9 11 13 15 17 19 21 23 25 27 29 
2*J 

Fig. 34. Contribution from If 7 proton pairs coupled to the maximum spin J = 6 

2 

(see text). 51 Mn minus 51 Fe. 

The sudden change of structure at J = y should also show up in the elec- 
tromagnetic transitions. Indeed, this is the case and the general trend of 
the transition values up to J = y breaks at J = y dropping to nearly 

zero. For the El transitions we have B(E2,^-~ ->■ y 1- ) = 0.013 e 2 fm 4 , 
B(E2,%~ -> f ") = 2.012 e 2 fm A and for the Ml transition B(M1,%~ -> 
f~) = 0.00002 fix, a reduction of several orders of magnitude that makes the 
state y isomeric. 

8 Magnetic moments of the N=28 isotones 

New measures of the magnetic moments of the N=28 isotones [43,55] have 
brought up again the question of the sensitivity of these observables to the 
effective interactions and to the valence space truncations as well as the rel- 
evance of the use of bare or effective (/-factors. We have examined this issue 
using the new interaction KB3G and we compare its results with those ob- 
tained using FPD6. 

The ground and low-lying states of the N=28 isotones are dominated by 1/t/ 2 
configurations. In the even-even cases, what is measured is the magnetic mo- 
ment of the first excited 2 + , therefore the closed shell nuclei 48 Ca and 56 Ni are 



40 



excluded from this systematics. In the I/7/2 limit, the value of the 2 + mag- 
netic moment, which is the same for 50 Ti, 52 Cr and 54 Fe, is 3.31 /j.^, using bare 
g-factors. The usual choice of effective g- factors, g^f—0.75 g^ are , g\=\.\ /in 
and g l u =-0.1 fiN, gives 3.08 hn- For the odd isotones, the ground states have 
J=7/2 and /i=5.79 /ijv with bare (/-factors and /i=5.39 /ijv with the effective 
^-factors. The differences amount to 7%. When mixing is allowed, the ten- 
dency is to reduce the difference between the results using bare and effective 
g-factors. 

In table 17 we have collected the results for the even nuclei calculated with 
KB3G and FPD6 at different truncation levels. For KB3G, a large reduction 
(20%) occurs between t=0 and t=3, while a marginal 5% more is obtained 
in the full calculation. For FPD6 the reductions are larger (30% and 10% 
respectively). The FPD6 results are too small compared with the experiment 
even if bare ^-factors are used, while those of KB3G are reasonably close to the 
measured values. The reason for the too large reduction that FPD6 produces 
can be related to the larger mixing with the H5/2 orbit, due to its very low 
location that we have already discussed. The experimental values of ji show a 
slight tendency to decrease with the number of protons, a trend that is weaker 
or even absent in the theoretical predictions. 
Table 17 

Magnetic moments of the first excited 2 + states of the even N=28 isotones (in ^1^). 





so Xi 


52 Cr 


54 Fe 


EXP 


2.89(15) 


2.41(13) 


2.10(12) 


KB3G (full) (bare) 


2.52 


2.50 




KB3G(t=3)(bare) 


2.65 


2.67 


2.56 


KB3G(t=3)(eff) 


2.49 


2.57 


2.52 


FPD6 (full) (bare) 


2.28 


1.90 




FPD6(fufl)(efT) 


2.12 


1.87 




FPD6(t=5)(bare) 


2.30 


2.08 


2.10 


FPD6(t=5)(eff) 


2.14 


2.04 


2.12 


FPD6(t=3)(bare) 


2.49 


2.38 


2.26 


FPD6(t=3)(eff) 


2.32 


2.32 


2.29 



In table 18 we present the results for the odd isotones. In this case KB3G 
seems to reproduce the decreasing trend better than FPD6, that gives a flatter 
behaviour. Considering the expected reductions in going from t=3 to the full 
calculation, the agreement can be considered quite decent, although in this 
case the use of effective ^-factors worsens it. As in the case of the even isotones, 
FPD6 gives too small values. 
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Table 18 

Magnetic moments of the 7/2~ ground states of the odd N=28 isobars (in (j,n)- 





49 Sc 


51y 


53 Mn 


55 Co 


EXP 


- 


5.15 


5.02(1) 


4.822(3) 


KB3G (full) (bare) 


5.34 


4.99 






KB3G(t=3)(bare) 


5.30 


5.05 


4.91 


4.75 


KB3G(t=3)(eff) 


5.00 


4.80 


4.71 


4.61 


FPD6(t=3)(bare) 


5.09 


4.89 


4.88 


4.81 


FPD6(t=3)(eff) 


4.77 


4.65 


4.67 


4.63 



9 Conclusions 



In this work we have extended our previous full p/-shell studies of the A=48 
and A=49 isobars three units of mass. In order to be able to treat properly 
the N=Z=28 shell closure and its surroundings we have introduced a mass 
dependence in the interaction KB3 and refined its original monopole changes. 
This results in the KB3G interaction. G emphasizes that the new interaction 
produces the right quasiparticle gaps for all the N=28 isotones. With this 
interaction -that can be interpreted as a natural extension of KB3, equivalent 
to it for the nuclei studied already- we achieve an extremely high quality 
description of the abundant experimental data available in the mass region 
A=50/51/52. Besides, we are now in a much better situation concerning our 
predictive power for higher masses. We have also studied the beta decay of the 
isotopes of Ca, Sc, V and Ti, computing half-lives and strength distributions, 
that agree nicely with the experimental information too. In the final sections 
we dwell on the magnetic moments of the N=28 isotones and on the Coulomb 
displacement energies along the yrast band of the A=51 mirror nuclei. 
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